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We present a first-principles atomistic thermodynamics framework to describe the structure, com- 
position and segregation profile of an alloy surface in contact with a (reactive) environment. The 
method is illustrated with the application to a Ag3Pd(lll) surface in an oxygen atmosphere, and 
we analyze trends in segregation, adsorption and surface free energies. We observe a wide range 
of oxygen adsorption energies on the various alloy surface configurations, including binding that 
is stronger than on a Pd(lll) surface and weaker than that on a Ag(lll) surface. This and the 
consideration of even small amounts of non-stoichiometries in the ordered bulk alloy are found to be 
crucial to accurately model the Pd surface segregation occurring in increasingly O-rich gas phases. 

PACS numbers: 68.43.Bc,71.15.Mb,81.65.Mq 
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I. INTRODUCTION 

The chemical properties of an alloy surface like its cat- 
alytic function or corrosion resistance depend intricately 
on the detailed surface composition. This compositional, 
and concomitant geometric and electronic structure is 
not only different from the corresponding surfaces of the 
parent metals, but differs often also substantially from 
the bulk alloy due to segregation of one species to the sur- 
face. The resulting segregation profile can depend on the 
surface orientation, the bulk composition, and tempera- 
ture. A reactive environment, such as one where cataly- 
sis or corrosion occurs, can further modify the chemical 
composition of the surface, if one alloy component in- 
teracts more strongly vifith a gas phase species than the 
others, resulting in adsorbate-induced segregation. It is 
consequently usually not possible to predict the prop- 
erties and functions of alloy surfaces in materials science 
applications solely on the basis of the formal bulk compo- 
sition; the effects of the corresponding (possibly reactive) 
environment must explicitly be considered. 

Previous theoretical attempts to address the segre- 
gation thermodynamics at an alloy surface from first- 
principles have mostly focused on the interdependence 
with the bulk alloy reservoir. This includes tabulated 
segregation energies of a transition metal impurity in 
another transition metal^'^, or approaches based on the 
coherent potential approximation"^, screened generalized 
perturbation method"* or surface cluster expansions'"^. 
The understanding emerging from these studies points 
at a crucial influence of the bulk reservoir structure and 
composition on the surface segregation profile, even if 
this concerns only slight excesses of one type of atom 
in samples with equal nominal bulk composition. This 
complexity has not yet been met in the scarce attempts 
to also explicitly account for the effect of a (reactive) 
environment^"*^, but the results reveal nevertheless an 
at least equally sensitive dependence of the adsorbate- 
induced segregation on the detailed adsorption properties 



and gas phase conditions. 

In view of these findings we set up a first-principles 
atomistic thermodynamics*'^"*^ framework based on 
density-functional theory that is able to address the 
structure and composition of an alloy surface while quan- 
titatively accounting for the dependencies on both the 
bulk alloy and the gas phase reservoir. The method is 
applied to obtain the segregation profile of a Ag3Pd(lll) 
surface in contact with an oxygen atmosphere, and we 
illustrate the intricate role of the bulk alloy reservoir by 
discussing the effect of already a small amount of non- 
stoichiometries (in form of anti-sites) in a nominally or- 
dered AgsPd bulk alloy. Not surprisingly we find Ag- 
segregation to occur under 0-poor environmental con- 
ditions, e.g. typical for ultra-high vacuum experiments, 
while Pd-segregation is preferred with increasing oxygen 
content in the gas phase. Disentangling the various fac- 
tors contributing to this segregation profile reveals, how- 
ever, an intriguing complexity that needs to be captured 
in the modeling to obtain quantitatively accurate results. 
In the course of this analysis we can thus scrutinize var- 
ious approximate treatments suggested in the literature, 
like the neglect of non-stoichiometries in the second or 
deeper substrate layers, the substitution of bulk segrega- 
tion with that of dilute impurities in the parent metal 
hosts, or the interpolation of bonding properties at the 
alloy surface from those known at the parent metal sur- 
faces. 



II. THEORY 

A. Surface free energy 

In the following we describe a first-principles atomistic 
thermodynamics framework to address alloy surface seg- 
regation in reactive environments. Our modeling aims at 
a system composed of a gas phase containing species j at 
partial pressures pj and a common temperature T, and 
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a multi-component alloy ^„^(b„,k)-5«B(,„,„ • • • /n.(„„,k) > 
where a::i(buik) = "i(buik)/ Z^j "i(buik) is the mole fraction 
of species i {i = A, B, /) in the alloy bulk. To evaluate 
the surface free energy of a particular surface of this al- 
loy in the given environment, we can concentrate on the 
part of the system that is affected by the surface, which 
contains (normalized to surface area Agurf ) Nj gas phase 
particles and a total of A'' solid particles at Xi mole frac- 
tions. The most stable surface striicturc and composition 
is then the one, which minimizes the surface free energy, 
defined as 

7(T,te},{iV,},7V,{xO) = 
G{T,{pj},{Nj},N,{xi}) - 

^{Nxi)m-^NjHj 



-^surf 



where G{T,pj,Nj,N, Xj) is the Gibbs free energy of the 
surface, and /i^ and iij are the chemical potentials of alloy 
species i and gas phase species j, respectively. 

For clarity, we will develop the methodological frame- 
work in this work for a binary alloy in contact with a 
single component gas phase, with obvious generalizations 
to more component systems. Specifically, we consider a 
^S«Ag(buik)P^"pd(buik) alloy and an oxygen atmosphere de- 
scribed by an oxygen pressure po^ ■ In view of the total 
energy calculations described below, we will furthermore 
consider the surface affected part of the alloy to be de- 
scribed by an inversion-symmetric slab geometry expos- 
ing two identical surfaces. With Agurf then denoting the 
area of the surface unit-cell, eq. (1) simplifies to 

7(^,P02,-^o,-/Vsiab,a;Ag,a;pd) = (1) 

Gslab (T", PO2 > -^O , -^slab , XAg , Xpd) - 
.^slab^JAg^Ag - ^slaba;pdMPd " -^OMO 
Gslab {T,P02, No , Nsiah , XAg , Xpd) - 
-^slaba;Ag(/^Ag - UPd) - NglabfiPd - -^OMO 



2^surf 



2A 



surf 



where A^siab refers now to the total number of metal 
atoms per surface unit-cell in the slab, A^o similarly to 
the number of O atoms in the surface geometry, and in 
the second step we have exploited that the mole fractions 
sum to unity by construction (xAg + xpd = !)■ 

Within our thermodynamic theory we assume the sur- 
face to be in equilibrium with the surrounding gas phase, 
and correspondingly the oxygen chemical potential is de- 
termined in the entire system by the gas phase reservoir, 
Ho = l/2M02(gas)(I";P02)- The alloy surface is further- 
more assumed to be in equilibrium with the underlying 
bulk alloy reservoir. This means that the chemical po- 
tentials /XAg and /xpd are not independent either, but are 



related by the bulk alloy chemical potential /ibuik, i-e. the 
Gibbs free energy per atom of the bulk alloy. For a bulk 
alloy of composition Ag„^^^^^j^j Pd„p^(^^,^ , the relation- 
ship between the chemical potentials is therefore 

Mbulk = a;Ag(bulk)(MAg - MPd) + MPd , (2) 

where we have again used the constraint on mole frac- 
tions. With this equation and introducing the differ- 
ence in chemical potentials of the two metal species, 
A/iAg-Pd = MAg — MPd, we can rewrite eq. (2) to 
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Gslab(^5P02 5-^05-^slab,i^^Ag,^Pd) " 



(3) 



-^slabMbulk — -^slab A/iAg-Pd {^Ag — Xa^ 

No ,rr ^ 

-^M02(gas) {J- ,P02) 



(bulk) J 



This equation has now an intuitive structure, with the 

first two terms giving the free energy difference between 
the surface structure and an equivalent number of bulk 
alloy atoms. The third term accounts for a possible dif- 
ference in the surface and bulk stoichiometrics and the 
fourth term for the additional number of O atoms present 
in the surface geometry. Finally, we note that throughout 
this text, we use the sign convention that a more nega- 
tive Gibbs free energy will indicate a more stable state of 
the system. In the case of a gas phase chemical potential 
this translates to /Uo2(gas) approaching —00 in the limit 
of an infinitely dilute gas, since adding a particle will 
then yield an infinite gain in entropy. As a consequence, 
7 > indicates the cost of creating the surface between 
the solid bulk phase and the homogeneous gas phase. 



B. Dependence on gas phase and alloy bulk 
reservoir 

The influence of the surrounding gas phase on the sur- 
face structure and composition enters eq. (3) through 
the term No/2 fJo2(eas){T,Po2)- I^i increasingly 0-rich 
environments, this term will favor surface structures con- 
taining an increasing number of oxygen atoms per sur- 
face area. In our approach we assess this effect by 
computing the surface free energy for a range of oxy- 
gen chemical potentials. At sufficiently low /tio2(gas) 
(0-poor conditions), only clean surface structures with- 
out oxygen atoms will be stable, which then forms a 
natural lower bound for the considered range of oxy- 
gen chemical potentials. As an upper bound, we use 
A^02(gas)(0 - rich) = teas), where E},°'/^^^^ is the total 
energy of an isolated O2 molecule (including the zero- 
point energy). This roughly corresponds to conditions 
where oxygen would start to condense at the sample at 
sufficiently low temperatures^^. 
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At the accuracy level of interest to our study, the oxy- 
gen gas phase is well described by ideal gas laws, which 
then allows to readily relate the oxygen chemical poten- 
tial to specific temperature and pressure conditions. For 
this we write^^ 

/^O, (gas) (T, po. ) = i?o .teas) + 2 A^o (T, po, ) (4) 

= <teas) + 2AMo(r,p°) + ksT In , 

where ks is the Boltzmann constant, and A/xq contains 
now all temperature and pressure dependent free energy 
contributions of the O2 molecule internal degrees of free- 
dom. Aij,o{T,po2) can be obtained can be obtained 
within the ideal gas approximation, or alternatively, one 
can use tabulated enthalpy and entropy values at stan- 
dard pressure = latm to determine Afio{T,p^)^^. 
For an oxygen gas phase, both approaches yield virtu- 
ally identical results in the range of temperatures and 
pressures of interest to our study^®'^''. 

Apart from the dependence on the gas phase chemi- 
cal potential, eq. (3) shows that the surface free energy 
depends also critically on the difference of the chemi- 
cal potentials of the metal atoms A/ZAg-Pd- This vari- 
able reflects the change in bulk chemical potential due 
to changes in a;Ag(buik) > as can be seen by evaluating the 
derivative of eq. (2) 



Ik 



dx 



Ag(bulk) 



MAg - MPd + 



(5) 



^Ag(bulk) 



+ (1 - a;Ag) 



Ag(bulk) 
= MAg - MPd = A/iAg-Pd 



dx 



Ag(bulk) 



where the term in brackets is equal to zero by the Gibbs- 
Duhem equation^^. The value of A/iAg-Pd will therefore 
depend sensitively on the bulk reservoir structure and 
composition, e.g. whether the bulk reservoir is ordered 
or random, or whether there is a slight excess of one type 
of atom. This will vary from sample to sample, and can 
have a significant impact on the surface energetics, even 
for samples with equal nominal bulk composition^^. 

In our application below, we will demonstrate this 
crucial dependence on the detailed bulk reservoir struc- 
ture and composition by considering the case of just a 
small amount of non-stoichiometries (in form of anti- 
sites) in the ordered AgaPd bulk alloy. Rather than fo- 
cusing on specific data for a specific sample, the aim is 
thereby to discuss the general effects that such bulk anti- 
sites can have on adsorption-induced segregation. In- 
stead of explicitly computing A/XAg-Pd e.g. by means of 
cluster expansion techniques®'^'^""^'^, we therefore pro- 
ceed by roughly estimating the range that A/j.Ag-pd can 
span, and discuss the results below for the obtained two 
bounds. For the estimate we use as limits the phase 
separation of Ag, when the bulk reservoir is rich in Ag 
(Ag-rich limit), and the phase separation of Pd, when 
the bulk reservoir is rich in Pd (Pd-rich limit). Phase 



separation of Ag will occur, when /iAg > A*Ag.fcc, where 
/XAg,fcc is the chemical potential of Ag atoms in the fee 
bulk structure of the parent metal. Equivalently, phase 
separation of Pd will occur, when /xpd > AtPd.fcc, where 
MPd.fcc is the chemical potential of Pd atoms in the fee 
bulk structure of the parent metal. Using eq. (2) and 
approximating /ibuik with its value for the perfectly sto- 
ichiometric alloy sample /Ubuik(stoich) i we therefore arrive 
at 

AMAg-Pd(Ag-rich) < A/iAg-Pd < A/XAg-Pd(Pd-rich) (6) 

with 

MAgjfcc ~ Mbulk(stoich) \ 



A/iAg-Pd(Ag-rich) 



1 - a;Ag(bulk) / 



(7) 



and 



. / Mbulk(stoich) " MPd.fcc \ . ^ 

A/XAg-Pd(Pd-rich) — Z I [°) 

\ 1 - a;pd(bulk) / 

as a first estimate for the range of A^Ag-Pd due to non- 
stoichiometries in the ordered bulk alloy reservoir. 



C. Solid phase Gibbs free energies 

In order to evaluate the surface free energy of a specific 
configuration with eqs. (3) and (6) we need to compute 
the sohd phase Gibbs free energies, Gsiab, A*buik(stoich), 
A^Ag,fcc and /ipd.fcc- For this it is useful to decompose 
them into several contributing terms, namely 



G 



^total _j_ pvih j.^conf. 



(9) 



where is the total energy (excluding the zero point 

energy), _F^'*' the vibrational free energy (including the 
zero point energy), and the configurational en- 

tropy. A crucial aspect that governs our analysis of all of 
these terms is that the quantity of interest to us, namely 
the surface free energy, does not depend on absolute 
Gibbs free energies. What enters into eq. (3) is a dif- 
ference of Gibbs free energies of slab, bulk alloy and bulk 
metals. This can allow for quite some degree of cancella- 
tion, e.g. due to similar free energy contributions in the 
three systems or due to similar errors in the computed 
total energies. 

We start our analysis by noting that the pV term is 
completely negligible for the transition metal alloys dis- 
cussed here^^. As for the vibrational free energy, a proper 
evaluation would require a systematic computation of all 
vibrational frequencies at the surface and in the bulk 
systems^^. Here, it can be particularly beneficial to real- 
ize that only the differences of the vibrational free energy 
contributions matter for eq. (3), i.e. the surface free en- 
ergy is only affected by the changes of the vibrational 
modes. In the application discussed below, we will ad- 
dress the effect of on-surface O adsorption on the segre- 
gation profile of a Ag3Pd(lll) surface. In this case, the 
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vibrational modes of Ag and Pd atoms in the bulk alloy 
and bulk metal structures will be rather similar, as will 
be their vibrational modes at the close-packed surface. 
From the similarity of the calculated stretch mode of oxy- 
gen adsorbed at low coverage at Ag(lll) and Pd(lll), 
SOmeV^'* and GOmeV^^ respectively, we furthermore do 
not expect significant variations in the vibrational free 
energies of alloy surface structures with different metal 
composition. What will therefore mostly contribute are 
the changes of the vibrational modes of oxygen in the 
gas phase and adsorbed at the surface. To estimate this, 
we use a simple Einstein model eqiiivalent to the one 
employed in Ref. 16 and compute the vibrational free 
energy contribution when the characteristic frequency is 
changed from ~ 200 meV (O2 stretch frequency in the 
gas phase) to ~ 55meV in the adsorbed state. Allowing 
for O coverages up to 1 ML and varying the characteristic 
frequency in the adsorbed state by ±50 %, already the to- 
tal contribution to the surface free energy stayed always 
below 5meV/A^ for temperatures ranging up to 600 K. 
What matters for the discussion of 0-induced segregation 
is, however, the difference of this vibrational surface free 
energy contribution for different surface configurations. 
For surface configurations with the same O coverage (dif- 
fering only in the Ag and Pd concentration in the sur- 
face layers) this difference will be completely negligible in 
view of the above cited similarity of the 0-metal stretch 
modes. For surface configurations differing in O cover- 
age, the difference in the surface vibrational free energy 
will be a bit larger, yet still negligible compared to the 
total energy terms discussed below which are of the order 
of tens of meV/A^. Correspondingly, we will neglect the 
zero point and vibrational free energy contribution to the 
surface free energy in this study, but emphasize that this 
is not a general result, and the effect of vibrations needs 
to be carefully assessed in each application. 

This leaves as remaining term the configurational en- 
tropy. A full evaluation of this contribution is compu- 
tationally involved, since it requires a proper sampling 
of the huge configuration space spanned by all possible 
bulk and surface structures. Although modern statisti- 
cal mechanics methods like Monte Carlo simulations^^'^^ 
are particularly designed to efficiently fulfill this purpose, 
they still require a prohibitively large number of total 
energy evaluations to be directly linked with electronic 
structure theories^^. A way to circumvent this problem 
is to map the real system onto a simpler, typically dis- 
cretized model system, the Hamiltonian of which is suffi- 
ciently fast to evaluate. Cluster expansions (also termed 
lattice-gas hamiltonian approach)®'^'^°"^^ are a promi- 
nent example for such a mapping approach, where the 
considered system would in practice e.g. be described 
by a lattice of possible on-surface adsorption sites for 
the gas phase species"'^'''^^. Here, we will instead con- 
centrate on screening a number of known (or possibly 
relevant or representative) ordered surface structures by 
directly comparing which of them turns out to be most 
stable under which (T,po2)"Conditions, i.e. which of 



them exhibits the lowest surface free energy. For suffi- 
ciently low temperatures, the remaining configurational 
entropy per surface area in Gsiab is then only due to a 
limited number of defects in these ordered surface struc- 
tures and can be estimated to be below 2 meV/A^ for any 
T < 600 K^^. We will restrict the discussion of the ap- 
plication to Ag3Pd(lll) below to the temperature range 
well below this limit, where then also the consideration 
of an ordered AgsPd bulk structure with a small amoimt 
of anti-sites is justified^". Focusing on this temperature 
range we will neglect the configurational entropy contri- 
bution to Gsiab, as well as the small configurational en- 
tropy of the ordered stoichiometric alloy and metal bulk 
reservoirs (entering Atbuik(stoich) . MAg,fcc and /xpd,fcc)- The 
only remaining temperature dependence in our model 
comes then from the entropy of the oxygen gas phase (en- 
tering Aiio{T,p)), as well as somehow implicitly through 
the considered bulk alloy chemical potential range, cf. eq. 
(6) . This will not affect the conclusions drawn below, but 
we note that the additional limitation that comes with 
the then resulting direct screening approach is that its 
reliability is restricted to the number of considered con- 
figurations, or in other words that only the stability of 
those structures plugged in can be compared. Its pre- 
dictive power extends therefore only to those structures 
that are directly considered, i.e. the existence of unan- 
ticipated surface geometries or stoichiometrics cannot be 
predicted. With this in mind, the here employed direct 
screening approach to first-principles atomistic thermo- 
dynamics can still be a particularly valuable tool, since 
it allows, for example, to rapidly compare the stability 
of different structural models without the need to have 
them mapped onto a common lattice. In this sense, we 
view the present approach as complementary to a cluster 
expansion study, since it provides first insight into the 
possible effects of the surrounding gas phase by screen- 
ing a wide range of possible structures. If more detailed 
insight is required (e.g. the exact surface concentration 
profile), this can then selectively be refined by a cluster 
expansion study, in particular when phase coexistence or 
order-disorder phenomena at elevated temperatures are 
of interest. 

In the present work, we therefore approximate the dif- 
ference of Gibbs free energies entering eq. (3) by the cor- 
responding difference of the leading total energy terms, 
and our final working equation reads 

l{T,po2, No, Nsiah,XAg,Xpd) ^ (10) 

— ^— £;*i°^b"'(^o,iVsiab,arAg,arpd) - 

■^^surf . 

-^slab-E'bu*k(stoich) - -^slabA/iAg-Pd(a;Ag " a;Ag(bulk)) " 

-^^^o*te-)-^oAMo(T,poJ . 

What is then required in practice to evaluate the surface 
energy of a given surface geometry of Ag3Pd(lll) in a 

given oxygen gas phase characterized by A/uq {T, po^ ) are 
total energy calculations of £^*[°*^'(A^o, -^siab) XAg,xpd) of 
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the surface (in the supercell geometry), of ^-bu'i^stoich) 
of the stoichiometric AgaPd alloy bulk, and of E^^^^^^s^ 
of an isolated O2 molecule. The considered range for 
the bulk alloy chemical potential difference A/it^gi^'p^ is 
given through eq. (6), for which (with the equivalent 
approximations) total energy calculations of the parent 
Ag and Pd bulk metals {E^^g'^L, -E^mJcc) are additionally 
required. If we define the formation energy for the bulk 
alloy as 

^^/(Ag„.,„,,,,Pdnp,,,.„,) = 

TTitotal ptotal 
-'^bulk(stoich) ~ 2;Ag(bulk)-C/Ag,fcc ~ 

(1 - a;Ag(bulk))^^Pd,fcc ; (11) 

then this range spanned by A/^Ag-Pd is given by 

^MAg-Pd(Ag-rich) — ^MAg-Pd(Pd~rich) ~ 

^^/(Ag„.,<,„,,Pd„p,,,,,,,) 




III 



IV 



2;Ag(bulk) (1 — a;Ag(bulk)) 



(12) 



D. Total energy calculations 

We compute all total energies within pseudopoten- 
tial plane wave density-functional theory (DFT)'^^, using 
the generalized gradient approximation (GGA)^^ for the 
exchange-correlation functional. The employed Pd ultra- 
soft pseudopotential (USPP)^'^"'^^ has a 4(f5s^ electronic 
configuration and is based on a Koelling-Harmon rela- 
tivistic all-electron calculation^®. The cutoff radii were 
set at 2.5, 2.5, 2.0 bohr for the s, p and d channels. The 
Ag USPP has a 4(i^"5s^ electronic configuration, with 
cutoff radii set at 2.5, 2.5, 2.1 bohr for the s, p and d 
channels. For both elements the / angular momentum 
channel is used as the local potential, with two s, p and 
d projectors, and a non-linear core correction^^ is used 
with rcoro = 1-2 bohr. The O USPP has a 2s^2p^ elec- 
tronic configuration, with cutoff radii set at 1.3 bohr for 
the s and p channels. The d angular momentum channel 
is used as the local potential, with two s and p projectors, 
and a non-linear core correction with rcorc = 0.5 bohr. 
We found it advantageous to use a less accurate oxy- 
gen USPP with a larger partial core radius (rcore = 0.7 
bohr) to pre-relax our geometries at a lower computa- 
tional cost. Substituting the more accurate O USPP in 
a final calculation, we found the maximum forces to be 
very close to or less than the force criteria set in the ini- 
tial optimization (< 0.05 eV/A). In several cases where 
the maximum forces exceeded the optimization criteria, 
we reoptimized the geometry, but found that the sur- 
face energies changed by less than ImeV/A^, which is 
accepted as a negligible error in the remainder of the cal- 
culations. A plane wave cutoff of 350 eV was employed 
in all calculations using the less accurate O USPP, while 
we employed the double grid method'^'* in the calculations 
using the more accurate O USPP. The higher energy cut- 
off (700 eV) with which the density is treated in the lat- 



FIG. 1: Top views of the atomic geometry of five (2 x 2) surface 
unit-cell arrangements considered in this study. The shown 
configurations have a bulk like second layer composition, but a 
first layer composition that varies from pure Ag (I) to pure Pd 
(V). White colored spheres correspond to Ag atoms, whereas 
grey colored spheres correspond to Pd atoms. The larger 
spheres correspond to the first layer, and the smaller ones 
correspond to the second layer. 



ter method increases the computational time moderately, 
but resulted in smooth convergence behavior"^^. 

In order to ensure maximum compatibility between 
the bulk metal and bulk alloy calculations, we employed 
the same bulk unit-cell to describe both the LI2 AgaPd 
and the fee Pd and fee Ag structures. (20 x 20 x 20) 
Monkhorst-Pack k-point grids were then used to op- 
timize the lattice constants to ao(Ag3Pd) = 4.09 A, 
ao(Pd) = 3.93 A, and ao(Ag) = 4.09 A. Using these op- 
timized lattice constants, we obtain a formation energy 
A£;/(Ag3Pd) as defined in eq. (11) of -0.21 eV per for- 
mula unit. The (111) alloy surfaces are modeled in a 
periodic supercell geometry, employing 7-layer inversion- 
symmetric slabs with two identical surfaces. The outer 
two layers on each side of the slab were fully relaxed, 
keeping the inner three layers fixed. A vacuum region 
of 14 A ensures the decoupling of consecutive slabs. The 
Brillouin zone integrations were performed with (8x8x1) 
Monkhorst-Pack k-point grids for the (2 x 2) surface unit- 
cells. For the oxygen adsorption calculations, oxygen was 
placed on each side of the slabs, maintaining inversion 
symmetry in all cases. For the computation of Eq^^^^^^ 
we employed a (12 x 13 x 14) bohr cell and obtain an O2 
binding energy of — 5.91eV. From test calculations with 
up to twice as many k-points and an increased plane wave 
cutoff of 500 eV we estimate the numerical accuracy of the 
reported surface energies to be within 2meV/A^, which 
does not affect any of the physical conclusions drawn. 



III. RESULTS 

A. O-induced Pd segregation at Ag3Pd(lll) 

We illustrate our approach by addressing the O- 
induced segregation profile at a Ag3Pd(lll) surface. At 
the nominal stoichiometry each (111) layer in the as- 
sumed LI2 bulk crystal structure'^*' contains 3 Ag atoms 
and 1 Pd atom in a (2 x 2) arrangement. We account 
for a possible segregation at the surface by considering 
(2 x 2) surface unit-cell atomic configurations in which 
the composition in the topmost two layers is systemati- 
cally varied from pure Ag to pure Pd. This leads to 25 
possible substrate configurations, five of which are shown 
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FIG. 2: (Color online) Surface free energy of Ag3Pd(lll) in 
equilibrium with a Pd-rich AgsPd bulk reservoir, cf. eq. (8), 
as a function of oxygen chemical potential. Each line corre- 
sponds to one of the tested surface configurations, and only 
the few configurations that result as most stable for a range 
of oxygen chemical potential are drawn as dark (red) lines. 
Additionally shown as insets are top views of the most stable 
surface configurations in the same form as in Fig. 1, with 
adsorbed O atoms shown as dark small circles, Ag atoms as 
white circles and Pd atoms as grey circles. The dependence 
on the oxygen chemical potential is translated into pressure 
scales using eq. (4) for T = 300K and T = 600K (upper 
a;- axes). 



in Fig. 1. In addition to these clean surface structures, 
oxygen adsorption was then examined by placing from 
one up to four O atoms per surface unit-cell into the on- 
surface fee sites. This leads to a pool of 473 inequivalent 
surface structures (including some larger unit cells), in 
which the oxygen coverage varies from to 1 monolayer 
(ML) in 1/4 ML steps in the (2 x 2) surface unit-cells. 

Evaluating eq. (10) for all configurations we arrive at 
the surface free energy plots shown in Figs. 2 and 3. 
Figure 2 corresponds to the limiting case of a Pd-rich 
AgsPd bulk reservoir and Fig. 3 to the opposite limit 
of a Ag-rich AgsPd bulk reservoir as discussed in Sec- 



FIG. 3: (Color online) Surface free energy Ag3Pd(lll) in 
equilibrium with a Ag-rich AgsPd bulk reservoir, cf. eq. (7), 
as a function of oxygen chemical potential. Each line corre- 
sponds to one of the tested surface configurations, and only 
the few configurations that result as most stable for a range 
of oxygen chemical potential are drawn as colored lines. Ad- 
ditionally shown as insets are top views of the most stable 
surface configurations in the same form as in Fig. 1, with 
adsorbed O atoms shown as dark small circles, Ag atoms as 
white circles and Pd atoms as grey circles. The dependence 
on the oxygen chemical potential is translated into pressure 
scales using eq. (4) for T = BOOK and T = 600K (upper 
a;-axes) . 



tion IIB, cf. eq. (6). Within our thermodynamic theory, 
at each oxygen chemical potential the surface configu- 
ration yielding the lowest surface free energy (i.e. the 
lowest line in the figures) results as "most stable". Fo- 
cusing first on the Pd-rich situation, we observe different 
stable configurations in different oxygen chemical poten- 
tial ranges. At the lowest chemical potentials considered, 
the lowest surface free energies are naturally obtained for 
the clean surface structures with varying chemical com- 
position in the topmost two layers. Among those, the 
"most stable" surface configuration under these 0-poor 
environmental conditions (i.e. the lowest horizontal line 
in Fig. 2) corresponds to a Ag-terminated alloy, with 
100% Ag in the first layer and 100% Pd in the second 
layer. As expected'^^, we therefore obtain segregation 
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of the more noble metal to the surface under gas phase 
conditions that at room temperature are representative 
of ultra-high vacuum experiments, cf. the pressure scales 
in Fig. 2. 

Due to the last term in eq. (10), increasingly O- 
rich surface structures get more favorable with increas- 
ing A^Q. The higher the O concentration (A^oMsurf) at 
the surface, the steeper the slope of the lines in Figs. 2 
and 3. In the case of the Pd-rich bulk limit, a first oxy- 
gen containing surface structure gets more stable than 
the clean Ag-segregated surface above A/xq = — 0.89eV. 
This configuration corresponds to an oxygen coverage of 
a 1/4 ML and contains now already 75% Pd atoms in the 
first layer. At even higher O chemical potentials further 
surface structures with higher oxygen coverages get stabi- 
lized, each time then containing a first layer composition 
of only Pd atoms. Induced by the oxygen adsorption 
is thus a segregation of the more reactive metal to the 
surface, which effectively reverses the segregation profile 
corresponding to room temperature ultra-high vacuum 
gas phase conditions. 

While this intuitive trend comes out clearly from our 
theory accounting for the effect of the environment, the 
exact surface concentrations as reflected by the obtained 
"most stable" structures have to be considered with care. 
Our approach considers a finite set of ordered surface 
structures and neglects the configurational entropy con- 
tribution to the solid phase Gibbs free energies as dis- 
cussed in Section IIC. The closeness of the lowest energy 
lines in Fig. 2 reflects a number of structures that are just 
a little bit less stable at a given A/xq than the highlighted 
"most stable" structure. This indicates already that a 
more exhaustive search of the huge configuration space 
e.g. by means of a cluster expansion technique would 
likely reveal new ground state structures with composi- 
tions that are intermediate to those considered here and 
that exhibit a certain degree of disorder at finite temper- 
atures. 

While a quantitative determination of the segregation 
profile is therefore outside the scope of the present work, 
the direct comparison of Figs. 2 and 3 nicely illustrates 
the entangled dependencies of the surface composition 
on both the bulk alloy and the gas phase reservoir. The 
results for the Ag-rich bulk alloy reservoir summarized 
in Fig. 3 also indicate an 0-induced Pd segregation at 
increasing A/io- One major difference to the afore dis- 
cussed Pd-rich bulk alloy reservoir is that the stabiliza- 
tion of these 0-containing Pd-enriched surface structures 
occurs at higher oxygen chemical potentials. This obvi- 
ously reflects the higher cost at which the additional Pd 
surface atoms have to be taken out of the Ag-rich alloy 
bulk reservoir, and correspondingly a higher driving force 
from the gas phase is required to establish the same Pd 
content in the surface layers as compared to the Pd-rich 
bulk alloy reservoir. A closer look at Figs. 2 and 3, and 
equivalently inspection of eq. (10), reveals, however, that 
the effect of the different bulk alloy reservoir is not just a 
mere shift of all surface free energy curves to higher A/^o 



by an amount that corresponds to the range of A/iAg-Pd 
between our two considered bulk alloy reservoir limits, 
cf. eq. (12). Within our computational setup this range 
amounts to 0.28 eV, and this is indeed roughly the order 
of magnitude by which the lines in Fig. 3 are shifted 
compared to those in Fig. 2. Still, due to the prefac- 
tor (xAg — a;Ag(buik)) of the A/XAg-Pd-term in eq. (10) 
surface energies of different surface configurations are af- 
fected differently by the change in the bulk reservoir. Ac- 
cordingly, we observe also a different sequence of "most 
stable" structures in Fig. 3 compared to Fig. 2, which 
underscores the need to explicitly account for both the 
effect of the environment and the bulk alloy structure to 
describe the structure and composition of alloy surfaces 
in realistic applications. 

Another intricacy that is nicely illustrated by our data 
is the effect of compositional changes in the deeper sub- 
strate layers. As apparent from Figs. 2 and 3 the first- 
layer Pd-segregation induced in more oxygen-rich envi- 
ronments is accompanied by substantial changes of the 
second layer composition. These changes have quite some 
impact on the stability of the various surface structures 
as is revealed when restricting our structure data base to 
those structures with nominal AgaPd second layer stoi- 
chiometry. When evaluating eq. (10) for this subset of 
surface structures we obtain as "most stable" clean sur- 
face for the case of a Pd-rich bulk reservoir a structure 
that contains only 75% Ag in the first layer (in contrast 
to the 100% Ag termination when also considering non- 
stoichiometries in the second substrate layer) . While the 
remaining sequence of "most stable" structures with in- 
creasing Pd content in the first layer is similar to the one 
shown in Figs. 2 and 3 respectively, the oxygen chemical 
potentials at which the transitions between the differ- 
ent structures occur are shifted by sometimes more than 
0.1 eV compared to the values exhibited in Figs. 2 and 
3. As indicated by the drawn pressure scales in the fig- 
ures this can correspond to several orders of magnitude in 
pressure in the temperature range discussed here, from 
which we conclude that a quantitative modeling of the 
Ag3Pd(lll) surface requires at least the consideration of 
non-stoichiometries in the topmost two substrate layers. 



B. Contributing factors to adsorbate-induced 
segregation 

The results presented in the previous section clearly 
show that the adsorption-induced segregation profile at 
the surface results from the compromise between the ten- 
dency to lower the surface free energy by segregating the 
more noble metal species to the surface, and the tendency 
to lower the surface free energy by achieving stronger ad- 
sorbate bonding at surfaces enriched with the more re- 
active metal species. In order to further qualify these 
two competing trends we define as a fundamental quan- 
tity describing the adsorbate bonding the average oxygen 
binding energy Eb{No,XAg,xpd) at a given surface con- 
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figuration containing Nq adsorbed O atoms per surface 
unit-cell and with a detailed surface composition in the 
surface layers characterized by mole fractions of a;Ag Ag 
atoms and xpd Pd atoms, 



2Nr 



.(13) 



As before this definition holds directly for symmetric 
slabs with adsorption on both sides, and with the cho- 
sen sign convention and zero reference in form of the gas 
phase O2 molecule an average binding energy E'b < 
indicates that dissociative adsorption is exothermic. 

Similarly we define as a quantity describing the sur- 
face segregation the segregation energy £^sog(a;Ag, 2:pd) as 
the cost to change the reference clean alloy surface with 
nominal bulk stoichiometric composition into an alter- 
native clean surface configuration with concentrations of 
XAg Ag atoms and xpd Pd atoms in the topmost layers 
(either by rearrangement of the surface atoms including 
atoms in sub-surface layers, or by segregation of atoms 
to or from the bulk reservoir), 



^scg(a;Ag, xpd) 



[£;*i°ab'(0,iV,iab,a;Ag,spd) 



-^slab'lO' A^slab, a:^Ag(bulk) j 2;pd(bulk) ) 

A^siabAi?Agl'pd(2^Ag - a:^Ag(buik))] 



(14) 



The factor 1/2 accounts as before for the fact that we are 
considering symmetric slabs where the segregation occurs 
equally at both surfaces, and within our sign convention 
a segregation energy ii^sog < indicates that segregation 
is exothermic. 

With these definitions, we can now recast our working 
equation for the surface free energy, eq. (10), into an 
intuitive structure. 



liT,po2iNo,Ns\i,h,XAg,xpd) = 
i?scg(a;Ag, a;pd) 



(15) 



7bulk-tcrm. + 



-^surf 

^;b(A^O,a^Ag,a^Pd) _ 

-^surf 



No 



■A^o(7',P02) 



where 7buik-term. is the surface free energy of the ref- 
erence clean alloy surface with nominal bulk stoichio- 
metric composition. This form of the equation shows 
clearly that the surface free energy can either be lowered 
by an energetically favorable segregation of metal atoms 
{Escg < 0), or by exothermic adsorption of iVo oxygen 
atoms (Eh < 0), where increasing oxygen chemical po- 
tentials in the gas phase further favor configurations with 
an increasing number of oxygen atoms at the surface. 
Having disentangled these various contributing factors, 
we now proceed to discuss our detailed results for each 
one separately. This enables us to carve out in detail the 
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FIG. 4: (Color online) Segregation energies per surface area, 
iSseg/^surf, for all clean surface configurations considered and 
as a function of the Ag concentration in the topmost layer. 
(Yellow) circles show the segregation energies for the Ag-rich 
bulk reservoir and (green) rhombes the segregation energies 
for the Pd-rich bulk reservoir. The different data points at 
the same topmost layer Ag concentration represent the vari- 
ations due to differing Ag concentration in the second layer, 
and the solid circles viz. rhombes correspond to the surface 
configurations with nominal bulk composition in the second 
layer. Additionally shown by the blue line is the segregation 
energy estimate obtained by considering the segregation of 
diluted Ag impurities in Pd and diluted Pd impurities in Ag 
(see text). 



complexity behind each of them, thereby bringing us into 
the position to scrutinize various approximate treatments 
that have been suggested in the literature. 



1. Segregation energy 

Figure 4 shows the segregation energies as defined in 
eq. (14) for the 25 clean Ag3Pd(lll) surface configu- 
rations computed. As expected, the segregation energy 
decreases overall largely with increasing Ag content in the 
topmost layer, which is the driving force behind the sta- 
bilization of the Ag-terniinated clean alloy surface struc- 
tures at the lowest oxygen chemical potentials in Figs. 2 
and 3. In detail, this variation of the segregation energy 
is, however, not simple, and depends much on the bulk 
alloy reservoir and the second layer composition. For the 
100% Ag-terminated surfaces, the lowest segregation en- 
ergy is for example obtained for a 0% Pd concentration 
in the second layer in the case of the Ag-rich bulk reser- 
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voir, but for a 100% Pd concentration in the second layer 
in the case of the Pd-rich bulk reservoir. 

This complexity needs to be captured in the modeling 
to obtain quantitatively accurate results for the segrega- 
tion profile, and thus dictates an explicit consideration 
of the structure and composition of both alloy bulk and 
surface. A prevalent alternative approach is instead to 
obtain estimates for the segregation energy by comput- 
ing the energy cost to segregate an impurity atom in a 
metal host to the surface^'^. In the present case, this 
would correspond to computing a diluted Ag impurity 
in the bulk or at the (111) surface of fee Pd, or a di- 
luted Pd impurity in the bulk or at the (111) surface 
of fee Ag. Within our computational setup, we obtain 
such impurity segregation energy estimates by calculat- 
ing the total energy difference when once placing a Ag 
atom in the central layer and once in the topmost layer 
of a seven layer (2 x 2) Pd(lll) slab (at the optimized Pd 
lattice constant), as well as respectively doing the same 
for a Pd atom and a Ag(lll) slab (at the optimized Ag 
lattice constant). The values obtained are — O.SleV for 
the case of a Ag atom in Pd, and -1-0.24 eV for the case 
of a Pd atom in Ag, which compare well to the equiva- 
lent numbers computed by Ruban and coworkers^ using 
the local-density approximation as exchange-correlation 
functional (-0.26 eV and -1-0.28 eV, respectively). 

Using these numbers the total segregation energy of 
a given surface configuration is then simply obtained by 
appropriately summing the impurity segregation costs. 
For a 100% Pd-terminated surface this corresponds e.g. 
to moving three Ag atoms per (2 x 2) surface unit-cell to 
the bulk and moving three Pd atoms per (2 x 2) sur- 
face unit-ceh from the bulk (-3 x (-0.31 eV) + 3 x 
0.24 eV = 1.65 eV). Normalized to the area of the (2 x 2) 
Ag3Pd(lll) surface unit-cell, this yields SsogMsurf — 
57meV/A^ in this simple model. The blue line in Fig. 
4 shows the corresponding data for such variations in 
the topmost layer Ag concentration, which reproduces 
the overall trend of the calculated DFT data surprisingly 
well (in particular when comparing to the filled symbols 
representing the DFT data where also the second layer 
composition has been kept at the bulk stoichiometry). 
By construction, the model can however not account for 
the dependence on the bulk alloy reservoir and leads thus 
to quantitative errors of the order of several lOmeV/A^ 
in the Egeg/Aguri contribution to the surface free energy. 
This is quite sizable considering the very magnitude of 
i?seg/^surf and 7 reflected in Figs. 2, 3 and 4, and under- 
scores the necessity for a more refined treatment of the 
segregation energy contribution which does account for 
the detailed alloy bulk and surface structure. 



2. Oxygen binding energy 

A similar complexity as for the segregation energy is 
observed, when we turn to the third term in eq. (15) 
and therewith to the term that is directly responsible for 
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first-layer Ag atoms 

FIG. 5: (Color online) Average binding energy per O atom 
as a function of the number of directly coordinated first-layer 
Ag atoms at a coverage of 9 = 0.25 ML, i.e. with one O atom 
per (2 X 2) surface unit-cell. For the considered adsorption 
into fee hollow sites at the close-packed (111) surface, this 
number of directly coordinated Ag atoms can range from zero 
to three. Additionally shown as (blue) squares is the average 
binding energy at this coverage at the fee sites of a Pd(lll) 
surface (corresponding to zero Ag atom coordination) and of 
a Ag(lll) surface (corresponding to three Ag atom coordi- 
nation). The (blue) line is a guide to the eye, representing 
a simple linear average between the binding energy at these 
two parent metal surfaces. 



the oxygen-induced segregation. Starting the discussion 
with the lowest O coverage that can be treated within 
our (2 X 2) surface unit-cells, i.e. 0.25 ML with 1 O atom 
per surface unit-cell. Fig. 5 shows the variation of the 
average binding energy with the number of first-layer Ag 
atoms to which the adsorbed O atom is directly coor- 
dinated. As expected, less exothermic binding energies 
are obtained with an increasing number of Ag atoms im- 
mediately involved in the oxygen bonding. This is in 
line with the known stronger bond strength of oxygen at 
the more reactive Pd(lll) surface compared to the more 
noble Ag(lll) surface. However, the situation at the 
alloy surface is not a mere linear composition weighted 
average between these two limiting cases as also shown 
in Fig. 5. In fact, we find average binding energies at 
Pd-rich Ag3Pd(lll) surface configurations that are even 
stronger than on Pd(lll), and average binding energies 
at Ag-rich Ag3Pd(lll) surface configurations that are 
even weaker than on Ag(lll). Judging from the different, 
known binding properties at the parent metal surfaces 
would therefore significantly underestimate the effect of 
oxygen-induced Pd surface segregation at this alloy sur- 
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face, i.e. also in this case a more simplified treatment that 
does not explicitly account for the alloy surface structure 
would not be appropriate. 

The range of average binding energies apparent in Fig. 
5 for the same number of directly coordinated first-layer 
Ag atoms is primarily due to the differing second layer 
compositions of the various surface configurations behind 
the data points. Interestingly, this effect of a varying 
number of Pd atoms in the second layer is different for 
different numbers of Pd atoms in the first layer. For 
the configurations with Pd-rich first layer compositions, 
increasing the number of Pd atoms in the second layer 
increases the average O binding energy, whereas for Ag- 
rich first layer compositions incorporation of Pd into the 
second layer substantially reduces the average O bind- 
ing energy. We find a similar effect of the second layer 
composition also at all other computed O coverages up 
to IML, each time leading to variations of Eh of the 
order of ~ 0.3 — 0.6 eV per O atom. For the most O- 
rich surface structures this amounts to variations of the 
term NoEh/Asuri in eq. (15) of up to ~ 80meV/A2, 
which is then of comparable order, if not larger than the 
corresponding variation of the i?scg/^surf term with non- 
stoichiometries in the second layer. The non-trivial in- 
terplay between these two equally-sized contributions to 
the surface free energy is thus what dictates at least the 
consideration of non-stoichiometries in the topmost two 
layers when aiming at a quantitative description of the 
adsorbate-induced segregation at this Ag3Pd(lll) sur- 
face. 

From our data, we also arrive at a first estimate of the 
impact that the third layer atoms have on the surface ad- 
sorption properties. For surface configurations with pure 
first and second layer compositions, i.e. either 100 % Ag 
or 100% Pd in either of the two layers, the four fee ad- 
sorption sites in our (2 x 2) surface unit-cells are com- 
pletely equivalent, except for variations in the type of 
atoms in the third bulk-stoichiometric layer (i.e. with 
three Ag atoms and one Pd atom per (2 x 2) cell). In 
these cases, we obtained variations of up to 0.1 eV per O 
atom in E^, depending on whether the oxygen atom was 
located directly over a Pd atom in the third layer or not. 
Compared to the average binding energy variations due 
to differing compositions in the topmost two layers, these 
are still noticeable, but already much smaller variations. 
On this basis, we would therefore conclude that a consid- 
eration of non-stoichiometries in the topmost two layers 
should capture the essential physics of oxygen-induced 
segregation at Ag3Pd(lll). 

Apart from the dependence on the surface composition 
the average O binding energy also displays a strong vari- 
ation as a function of coverage. Figure 6 illustrates this 
for the subset of computed configurations with either a 
100% Ag or a 100% Pd concentration in the topmost 
layer. In both cases, we observe a strong reduction in 
bond strength with increasing number of oxygen adsor- 
bates at the surface. For the Ag-terminated surfaces, 
this has the consequence that dissociative adsorption be- 
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FIG. 6: (Color online) Average binding energy i?b per O atom 
as a function of the oxygen coverage, ranging from 0.25 ML 
to 1 ML (corresponding to one to four O atoms per (2 x 2) 
surface unit-cell). Shown is only the data for surface config- 
urations with a 100 % Ag concentration in the topmost layer 
(yellow squares) or with a 100 % Pd concentration in the top- 
most layer (red circles). The different data points at the same 
coverage correspond to surface configurations with varying 
second layer composition. The full Ag second layer gives the 
lower energies and the full Pd second layer the higher energies. 



comes endothermic at all configurations with coverages 
6 > 0.25 ML, which is very similar to the findings at the 
pure Ag(lll) surface^". Within the scatter in the data 
shown in Fig. 6, which is due to the varying second layer 
composition of the surface configurations displayed, the 
decrease of the average binding energy with coverage is 
roughly linear for both Ag and Pd terminated surfaces. 
Again, this is similar to the bonding at the parent metal 
surfaces Ag(lll) and Pd(lll)^°. The amount, by which 
i?b decreases with coverage is, however, not the same for 
the two differently terminated surfaces. Whereas for the 
Pd-terminated surfaces the average binding energy de- 
creases by about 1 eV in the computed coverage range 
0.25 < 9 < IML, it decreases for the Ag-terminated 
surfaces by about 1.5 eV. In other words, the anyway 
stronger bonding at Pd-enriched surfaces at low O cover- 
ages becomes even more favored at increasing O coverage 
in comparison to Ag-enriched surfaces with the same oxy- 
gen load. This will further enhance the tendency towards 
Pd surface segregation at increasing O chemical poten- 
tials, and shows that a restricted information about low 
coverage bonding properties will also not be sufficient to 
correctly model the segregation profile in a reactive envi- 
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ronment like the presently studied oxygen atmospheres. 

3. Gas phase chemical potential 

Even though the average binding energies per O atom 
decrease with increasing oxygen coverage, the total con- 
tribution of A^o times the smaller in eq. (15) may 
still lead to lowest surface free energies for surface con- 
figurations with higher oxygen load. At increasing oxy- 
gen chemical potentials such structures are further fa- 
vored by the last term in cq. (15), i.e. the increasing 
abundance of oxygen in the gas phase helps to stabi- 
lize 0-rich surface structures. Within our calculated set 
of surface configurations, this leads to the sequence of 
"most stable" structures in Figs. 2 and 3, with increas- 
ing concentrations of adsorbed O atoms at higher values 
of A/io- At sufficiently high oxygen chemical potentials 
also more 0-rich surface configurations beyond the sim- 
ple on-surface adlayers considered in the present study 
may get stabilized. Apart from ultra-thin so-called sur- 
face oxide films, which have been characterized for both 
parent metal surfaces^^'^"'* , this will ultimately be thicker 
bulk-like oxide films, and due to their higher stability 
most likely palladium oxides. 

Assuming that what forms at the surface is simply the 
most stable palladium bulk oxide, PdO, we can estimate 
the oxygen chemical potential at which such a bulk oxi- 
dation occurs from the stability condition^^ 

MPdO < MPd + MO • (16) 

Within the same approximations as applied in Section 
lie, we then arrive in the limit of the Pd-rich bulk reser- 
voir (/ipd = Mpd(motal)) at 

^^bulkoxid.,Pd-rich ^ ^EfiPdO) , (17) 

where the formation energy of bulk PdO is defined as 

Ai?/(PdO) = i?^,°d*S'.bulk - ^^Pdlcc - l/2£;S°*feas) • 

(18) 

Similarly, one obtains in the limit of the Ag-rich bulk 
reservoir (/XAg = AiAg(motai)) 

^^buikoxid.,Ag-rich ^ A£;/(PdO) - A£;/(Ag3Pd) , 

(19) 

with the alloy formation energy as defined in cq. (11). 

Within our computational setup, we determine a value 
of Ef (PdO) = — 0.86eV per formula unit, so that the 
chemical potentials corresponding to the onset of bulk 
oxide formation are ^^^"i^°^"i->P'i-"<=h _ _q gg^y 
^^buikoxid.,Ag-rich ^ _o.65eV. In both cases, the onset 
of bulk oxide formation occurs therefore just at slightly 
higher oxygen chemical potentials than the first stabi- 
lization of surface configurations with adsorbed O atoms 
within the set of structures compared in our study (at 
A/xo = -0.89 eV and -0.68 eV respectively, cf. Figs. 2 



and 3). This is largely different compared to the situa- 
tion at Pd(lll), where a p(2 x 2) structure with 0.25 ML 
on-surface O adatoms in fee hollow sites gets stabilized 
at oxygen chemical potentials that are significantly lower 
than the onset of bulk oxide formation^'*. The reason 
for this difference is the high cost of segregating the 
scarce Pd in AgsPd to the surface that is not necessary 
at Pd(lll). The average O binding energy we compute 
at 0.25 ML coverage at the Ag-rich surface termination 
(100% Ag in the first layer, 100% Pd in the second layer) 
that is most stable at the lowest A/uq is only — 0.13eV/O 
atom, which should be compared to our computed av- 
erage binding energy of — 1.29eV/0 atom at the pure 
Pd(lll) surface at this coverage. Only by segregating 
Pd atoms to the first layer (and thus involving them 
into the direct coordination to the O adsorbate), does 
the Eh of corresponding Ag3Pd(lll) configurations come 
into the same range as the Pd(lll) binding energy, e.g. 
— 1.29eV/0 atom for the configuration with 100% Pd in 
the first layer and 100% Ag in the second layer, cf. Fig. 
5. However, stabilizing such configurations then involves 
the cost of segregating Pd to the topmost layer and thus 
delays their stability range up to oxygen chemical poten- 
tials, where bulk oxide formation is already about to set 
in. 



Within the set of calculated (2 x 2) surface unit-cells, 
we can obviously not exclude that there are not more di- 
lute Ag3Pd(lll) adsorption structures with 9 < 0.25 ML 
that could already become most stable at a A/io that is 
noticeably lower than the limit of bulk oxidation. How- 
ever, the essential physics should already be captured 
within our considered subset of configurations, and this 
is that Pd and Ag are chemically simply too similar, or 
more precisely that the 0-Pd binding is not stronger by a 
sufficient amount compared to the 0-Ag binding, to over- 
come the high Pd segregation cost in AgsPd and induce a 
noticeable adsorbate-induced segregation at already low- 
est oxygen chemical potentials. Hence, segregation of Pd 
only sets in, when the driving force from the gas phase is 
already so high, that this directly initiates the formation 
of bulk-like oxide films at the surface. In this respect, it 
will be interesting to investigate how this compares to the 
situation in more Pd-rich Pd-Ag alloys, or to the case of 
binary alloys formed of more chemically different metals, 
like Pt-Ru alloys. There the strong 0-Ru binding will 
exert a large driving force for Ru segregation to the top- 
most layer^°, which on the other hand needs to overcome 
the much higher segregation energy to enrich the surface 
with the more reactive Ru metal^'^. At which oxygen 
chemical potentials this then leads to the stabilization 
of Ru-rich 0-containing surface structures, and how this 
A/zo then compares to the onset of bulk oxidation, will 
depend sensitively on the availability of Ru atoms, and 
therewith on the detailed structure and composition of 
the bulk alloy reservoir. 
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IV. SUMMARY 

In conclusion, we have presented a first-principles 
atomistic thermodynamics approach to describe the 
structure and composition of an alloy surface in con- 
tact with a (reactive) environment. Accounting for both 
the dependency on the bulk alloy and the gas phase 
reservoir, this approach provides the appropriate frame- 
work for a detailed discussion of the driving factors be- 
hind adsorbate-induced surface segregation. We illus- 
trate this for the Ag3Pd(lll) surface exposed to an oxy- 
gen atmosphere, where we obtain an inversion of the 
segregation profile in increasingly oxygen-rich environ- 
ments. Whereas a minimal segregation energy stabilizes 
Ag-tcrminated surface structures in ultra-high vacuum 
type gas phase conditions, the much stronger oxygen 
bonding favors increasingly Pd-rich terminations in at- 
mospheres with higher oxygen content. Our analysis 
shows that the details of this transition are intricately 
coupled to the bulk reservoir structure, and already small 
non-stoichiometries in the nominally ordered AgsPd bulk 
alloy can not only change the oxygen pressures required 
for the transition by several orders of magnitude, but 
also afi'ect the explicit segregation profile itself. This 
highlights the restricted relevance of detailed segrega- 
tion data from ultra-high vacuum experiments for ap- 
plications in realistic environments, as well as the com- 
plexity that needs to be captured in the modeling to 
obtain quantitatively accurate results. In this respect 
our study shows that at least for the Ag3Pd(lll) sur- 
face simplified treatments based on impurity segregation 



and adsorption data of the parent metals are not able 
to correctly describe the adsorbate-induced segregation, 
nor that restricted information about the low-coverage 
adsorbate bonding at the alloy surface or the exclusive 
consideration of non-stoichiometries in the topmost sub- 
strate layer would be sufficient. This holds already for 
the case where the reactive environment leads only to 
the formation of on-surface adsorbate layers, and is even 
more pronounced when considering the possibility of a 
more complex restructuring like here the oxide forma- 
tion in an 0-containing environment. Also in the latter 
case the ruling factors are the availability of the more 
reactive metal species in the alloy bulk and the increased 
oxygen bond strength that this species can provide. In 
the studied case of the Ag3Pd(lll) surface the high cost 
of segregating the scarce Pd to the surface delays the 
stabilization of Pd-enriched surface structures in fact al- 
most up to gas phase conditions where the formation of 
thicker palladium oxide films would already become ther- 
modynamically favorable. We speculate that this will be 
different in more Pd-rich Ag-Pd alloys, or in alloys com- 
posed of species with a larger difference in their O binding 
capabilities. 
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